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Essential Electrodynamics: Solutions 


Preface 


Preface 


This book gives the solutions to the exercises at the end of each chapter of my book “Essential 
Electrodynamics” (also published by Ventus Publishing ApS). I recommend that you attempt 
a particular exercise after reading the relevant chapter, and before looking at the solutions 
published here. 

Often there is more than one way to solve a problem, and obviously one should use any valid 
method that gets the result with the least effort. Usually this means looking for symmetry in the 
problem - for example from the information given can we say that from symmetry arguments 
the field we need to derive can only be pointing in a certain direction. If so, we only need to 
calculate the component of the field in that direction, or we may be able to use Gauss’ law or 
Ampere’s law to enable us to write down the result. In some of these exercise solutions the 
simplest route to the solution is deliberately not taken in order to illustrate other methods of 
solving a problem, but in these cases the simpler method is pointed out. 

The solutions to the exercise problems for each chapter of “Essential Electrodynamics” are 
presented here in the corresponding chapters of “Essential Electrodynamics - Solutions”. 

I hope you find these exercises useful. If you find typos or errors I would appreciate you 
letting me know. Suggestions for improvement are also welcome - please email them to me at 
protheroe.essentialphysics@gmail.com. 

Raymond J. Protheroe 

School of Chemistry & Physics, The University of Adelaide, Australia 


Adelaide, May 2013 
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Electrodynamics and conservation laws 


1 Electrodynamics and conservation laws 


1-1 A magnetic dipole of moment m = mz is located at the origin. A thin circular conducting 
ring of radius a vibrates such that the position of its centre is r = [zo + bcos(cjt)]z with 
b « a « zq. The plane of the ring remains parallel to the x—y plane during the vibration. 
Find the emf around the ring in the </> direction. 

Solution 

The magnetic field of the dipole is 


B(r) 


y o 3r(m • r) — r 2 m 
47 t r 5 


( 1 . 1 ) 


Since 6 < a C zo we can approximate the magnetic field anywhere on the ring as it vibrates 
by 


B[z r (t)z)] 


MO 

47T 


3z r z (mz 


yo mz 
2i r z^ 5 


where z r (t) = [zo + &cos(c<;t)] is the height of the ring. 
The magnetic flux through the loop is 


= 7TQ? - [zq + bCQs(ujt)] 3 , 


2n 

2 


y^ma 

24 

yoma 2 

~^r 


-3 


14-cos (cut) 

z o 

b 

1 — 3— cos (uot) 

z o 


since a « Hence, 


£ = - 


d$B 

~dT' 

3yoma 2 buj 

2zq 


sin (uot). 


( 1 . 2 ) 


(1.3) 

(1.4) 

(1.5) 


( 1 . 6 ) 

(1.7) 


1-2 A thin disc of radius a and height h contains charge +q uniformly distributed throughout 
the disc. The disc is located with its centre at the origin, and rotates about the z-axis with 
angular velocity a? = uoz. 
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Electrodynamics and conservation laws 


(a) Using cylindrical coordinates but with R being the cylindrical radius to avoid confusion 
with the charge density p( r), specify the current density J(i?, </>, z) as a function of position. 
In the limit h <C a find the magnetic dipole moment. 

(b) Consider a circular loop of radius Rq around the z-axis at height zo above the disc for 
the case Rq a zq. Find the magnetic flux through the loop, and hence find the vector 
potential at the loop. 

(c) If, due to friction in the axle, the disc’s angular velocity is decreasing exponentially 

with time t as u{t) — where to is the decay time scale, find the electric field at 

the loop at time t — 0. 

Solution 


(a) Within the disc, pi r) = —and v(r) = Rucf), and so 

7 va z h 


J(r) = p(r) v(r) = p(r)Ru;4> = 


(1.8) 


The dipole moment is 


m 


= I/r x J(r) d 3 r, 

; ^ r 

a 2 Jo 


( qRoo 
na 2 h 


<fi ) 27 vRdR, 


R? dRz , 


qua 


z. 


(1.9) 

( 1 . 10 ) 

( 1 . 11 ) 

( 1 . 12 ) 


(b) The circular loop is close to the axis of the dipole, but a distance zo a away. The 
magnetic field of a dipole is 


B(r) = i^3[ 3 ( m -r)?- m ] 


B(0, 0, zq) -o 2m z 

4?rz o 


8ttZq 


qua z. 


(1.13) 

(1.14) 


The magnetic flux through the loop is then approximately 


= ttR 2 B(0,0,z 0 ) = qua 2 . 

O Zr\ 


(1.15) 
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We can obtain the vector potential from the magnetic flux using 


A • dr / B • dS — 

r Js 


(1.16) 


where loop V is the circular loop and S is any surface bounded by the loop. From symmetry 
arguments the vector potential must be in the <fi direction 


a/d \ 2 Mo a 2 Ro , 

A(R 0 ,(f),z 0 ) = = T^z—5-qu<p. 


2ttRq 


16n Zq 


(1.17) 


(c) If u>(t ) = ui^e tRo then 


F,(R 0 ,(j),z 0 ,t) = 


d 

-gj;M R o,<P,zo,t), 

Mo a2 Ro d -t/t 0 2 

~TZ -3 -<l u 0-EL e 

lo7r Zq at 

jio_a^RoquJo -t/to-2 
Kirr to ' 


(1.18) 

(1.19) 

( 1 . 20 ) 
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1-3 A light rigid rectangular circuit with resistance R has mass m attached to the middle 
of the lower side (width s), and the top side is suspended horizontally using frictionless 
bearings to form a simple pendulum of length h as shown in the diagram below. In the 
absence of a magnetic field the position of the pendulum mass would be described by 
r m (t) ~ h0Q cos(u;t) x where u = yfgjh. A uniform magnetic field B points in the verti¬ 
cally upward direction. 



(a) Assuming the position of the pendulum mass is still described by r m (t) ^ h 6 0 cos (cot) x, 
what is the magnetic flux <&#(£) through the circuit, and hence the emf as a function of 
time? Take the direction around the circuit indicated by the arrow to correspond to positive 
emfs and currents. 

(b) What is the force on the lower side of the circuit due to the magnetic field? What 
is the instantaneous work done by the pendulum against this force? Compare this with 
instantaneous power dissipated in the circuit? What are the consequences of the presence 
of the magnetic field for the motion of the pendulum? 

Solution 

(a) The magnetic flux through the circuit is 

$#(£) = x-r m (£)s£> = h0QCOs(ujt)sB (1.21) 


an so the emf is 


S{t) 


d&B 

dt 


uohsB0 0 sin (wt). 


The current is 


I{t) = £(t)/R = ((jjhsB0o/R) sin(ut). 


( 1 . 22 ) 


(1.23) 
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(b) The Lorentz force on the lower side of the circuit is 
F (t) = I f dr' x B, 

J lower side 

= (ujhsBOo/R) sin(cdt) (sy) x (B z), 

= (ujhs 2 B 2 9o/R ) sin (ut) x, 


(1.24) 

(1.25) 

(1.26) 


and since v(t) = dr/dt = — cjoHOq sin (out) x, we note that F (£) is in the opposite direction to 
v(t). The instantaneous work done by the pendulum against this force is 


-fmech(^) = -F (t)-v(t), 

= — [(cjhs 2 B 2 0Q /R) sin(cdt) x] • [— uoh 6 q sin(cdt) x], 
= (<jj 2 h 2 s 2 B 2 6 o/ i?) sin 2 (o;t) . 


Note that The instantaneous power dissipated as heat in the resistor is 


£ 2 

Pheatyt) ~ — 


[cd/is£>0o sin(cdt)] z 


(1.27) 

(1.28) 
(1.29) 


(1.30) 


consistent with conversion of the mechanical energy of the pendulum to heat. 

(c) The pendulum’s amplitude 6q will decay over time. To determine the rate of decay, we 
can consider first consider the average rate of energy loss of the pendulum which is equal 
to the time-average of the power dissipated as heat 


(-^heat) 


(uhsBOo) 2 

2R 


(1.31) 


and compare this with the total energy of the pendulum which is equal to its kinetic energy 
at x = 0 


Wtot = i^max = ^m(uh6o) 2 . 


(1.32) 


The energy will decay exponentially with timescale 


_ Wtot _ 1 _ ( \2 

“ </w> “ 2 miu,Mo) (uhsBetf 


mR 

^B2’ 


(1.33) 
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that is, W ta t{t) oc e Cinergy. Since the energy is proportional to 0 q the pendulum amplitude 
9q will decay exponentially as 9o(t) oc where tq = 2r energ y. 


1-4 Consider the section of a two-wire transmission line shown below. Show that the self¬ 
inductance per unit length for the case where D a is given by 


L = 


— In — 
7r a 


(1.34) 



Solution 


We need to calculate the magnetic flux between the two wires for a section of length b due 
to current / flowing in both wires. For just one wire the magnetic field is 
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B=^ 

2np 

and integrating between the conductors we get the flux for length b 


Mo/ f D ~ a dp 

2?r J a p 


^1 /D-a\ 

27 r \ a J 


(1.35) 


(1.36) 


Multiplying by 2 (to include the flux from both wires) and dividing by b and the current 
we get the inductance per unit length 


L 


— In 
7r 


D — a 



(D > a). 


(1.37) 


1-5 Consider a coaxial cable as an infinite cylindrical inductor and find the inductance per unit 
length. 

Solution 



Consider length a with current / in the z direction on the inner conductor and / in the 
—z direction on the outer conductor. To find the magnetic field, use Ampere’s law for a 
coaxial circular loop of radius p, for p\ < p < p 2 . 


2typ B(p,(j), z) = p/. 


.-. B{p,<p,z) 



(pi<p< pi). 


(1.38) 


The magnetic flux through the rectangular loop T of length a and width (p2 — Pi) spanning 
radii pi < p < p2 as shown is 





p/ 
27rp 


a dp 


pal 
2i r 


In 



(1.39) 
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Hence, the inductance per unit length is 

L = = — In (—) . 

la 2n \pi J 


(1.40) 


1-6 The diagram shows a parallel plate capacitor. Find the current /(t), and the displacement 
current density J d between moving capacitor plates, and check that the total displacement 
current is Id(I) = I(t). Neglect fringing effects. 



Solution 

From the capacitance of a parallel plate capacitor we get the charge on the positive plate 


Q( ( )=C«)V = ^. 


So there is a current flowing 

_ dQ _ £qAV dx _ £qAV 
dt [x(t)] 2 dt [x{t)] 2V 

The electric field between the plates is 


(1.41) 


(1.42) 


m 


V ^ 
——x. 

x(t) 


(1.43) 


The displacement current density between the plates is then 

(i- 

Multiply Jd by the plate area Az to get the total displacement current 

Id = = —£o -r- 

[x 


m 


;Av(t) 


(1.45) 


<9E 

J d = ^o^rr = 
dt 


V dx ^ 

~ g 0 r / i\i9 ~TT X = — ^0' 


V 


c(t)] 2 dt 


x{t)\ 
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which is the same as the current / in the wire. 


1-7 Consider a straight piece of wire radius a and length Az, along which current / is flowing. 
The potential difference between the ends is AV. Find the Poynting vector at the surface of 
the wire and use it to determine the rate at which energy flows into the wire, and compare 
the result with Joules’s law. 

solution 



A z 


When a current flows along a resistive wire, work is done which shows up as Joule heating. 
The energy flows from the electromagnetic field into the wire, and the rate at which energy 
is flowing can be determined from the Poynting vector. If the wire has potential difference 
AV between its ends, the electric field is parallel to the wire and given by 
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Ay 

E = —i (1,46) 

From Ampere’s law the magnetic intensity of current / is circumferential, and at the surface 
has the value 


H =2^- 


Then, the Poynting vector is 


S = 



A VI ^ 
2na A 


(1.47) 


(1.48) 


The Poynting vector points radially inward and so the energy per unit time passing in 
through the surface of length A z of the wire, which has area 2 tt a A z, is 



S -dS 


AVI 

2n a A z 


2n a A z = AVI, 


(1.49) 


where dT, is used here for the surface element to avoid confusion with the Poynting vector 
S. The current density has magnitude J — I/na 2 , so the Joule heating power of the wire, 
which has volume 7ra 2 Az, is 



E • J d 3 r 



(1.50) 


The current and voltage are constant, and so the magnetic and electric fields also are 
constant, 


d 

dt 




d 3 r = 0. 


(1.51) 


The result, Poynting flux into the wire (Eq. 1.50) equals Joule heating power (Eq. 1.49) 
plus rate of increase of held energy inside wire (Eq. 1.51), is consistent with Poynting’s 
theorem. 


1-8 A long solenoid carrying a time-dependent current / (t) is wound on a hollow cylinder whose 
axis of symmetry is the z-axis. The solenoid’s radius is a, and it has n turns per metre. 
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(a) Write down the magnetic intensity H(r, t) and magnetic field B(r, t) everywhere. What 
is the energy density in the magnetic field inside the solenoid? 

(b) Find the electric field E(r,£) everywhere using Faraday’s law in integral form. 

(c) Find the magnetic vector potential A(r,£) everywhere. 

(d) Find the Poynting vector S(r, £) inside the cylinder, and hence the energy flux into 
a section of the cylinder of length h and the rate of increase of energy density inside the 
cylinder. Compare this with the rate of increase of magnetic field energy inside length h 
of the cylinder. 

Solution 

(a) From earlier work on electromagnetism we should know that 


B (p, 4>, z) 


IMynJ z (p < a) 

0 (p> a) ’ 


H(p,0,z) 


nlz, (p < a ) 

0 (p > a) 

(1.52) 


The magnetic energy density inside the solenoid is 

u _ -EP_ _ (jMynff _ /x 0 (n/) 2 

2/i o 2/xo 2 


(1.53) 


(b) As the current increases, so does the magnetic field which is uniform and remains in the 
z direction. Thus, the induced electric field vector must be perpendicular to the z direction, 
and from symmetry arguments must be in the ±0 direction, i.e. E(r, t) = Hence, 

using Faraday’s law 


y E • dr — 

-it* 1 ” 

(1.54) 

2irpE<p(p,t) = 

j —(7rp 2 )pondI/dt. (p < a) 

1 —(7va 2 )pQndI/dt. (p > a) 

(1.55) 

E(r ,t) = 

j —(p/ion/2) dl/dt<f>. (p < a) 

) —(a 2 p J on/2p)dI/dt<fi. (p > a) 

(1.56) 



(1.57) 


(c) To find the vector potential we can use 


y*A(r, t)-dr = f s B(r, t) • dS. 
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From symmetry arguments, the vector potential must be in <p direction, and therefore we 
should use a concentric loop T of radius p giving 


271 -pA# 

A 


up 2 pqtiI (p < a) 
7ia 2 ponI (p > a) 

(ponIp/2) </> (p<a) 

(ponla 2 /2p) <fi (p > a) 


(d) The Poynting vector at p — a is 




E(a, </>, z) x H(a, </>, z), 


- ) x (nl z), 


a o dl ^ 

-2 mnI dt P ' 

a d , % 9 ^ 

P, 


(1.58) 

(1.59) 


(1.60) 

(1.61) 

(1.62) 

(1.63) 
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and is directed radially inwards. 

To find the energy flux going into a section of the solenoid of length h we multiply by the 
its surface area 2irah 



(1.64) 


which is just the rate of increase of magnetic energy inside length h of the solenoid. 


1-9 Using the Maxwell stress tensor find the pressure exerted on a perfectly absorbing screen 
by an electromagnetic plane wave at normal incidence. 

Solution 

Assume the screen is in the x-y plane and the EM wave is 



(1.65) 


Momentum conservation in electrodynamics is expressed by 



( 1 . 66 ) 


where n is the outward normal unit vector at surface S which bounds volume V, and 



(1.67) 


The volume V corresponds to the space defined by z > 0, so n = —z and e j • n ^ 0 for 
j z. Hence we need only calculate T xz , T yz and T zz . 

For the EM wave above T xz = T yz = 0 and 

T zz — — -/^qE-EH - B • B^j — —u (1.68) 

where u is the energy density of the EM wave. 

The pressure is given by the rate of change of the momentum of particles and fields inside 
V, so we need to integrate over a portion of the surface, i.e. some area A of the x—y 
plane, and with V bounded by A and extending in the +z direction. Then the force on 
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the contents of V is 


F = 



dS, 


= z [T zz z ■ (-z)] A 
= —z(—u)A 
= uAz. 


(1.69) 

(1.70) 

(1.71) 

(1.72) 


Hence, the pressure p — F / A — u is equal to the energy density in the EM wave. 

If we had not been asked explicitly to use the Maxwell stress tensor we could have arrived 
at the same result in a simpler way by noting that the momentum density is g = S/c 2 , 
and that for an EM wave S = uc\a. Pressure will be the magnitude of the momentum flux 
gc — (S /c 2 )c = u. 



t) 
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2 Electromagnetic waves in empty space and linear dielectrics 


2-1 Prove that the spherical waves given by 

( 2 . 1 ) 

are solutions of the 3D wave equation. 

Solution 

We need to show that this is a solution of the three-dimensional wave equation 


f(*,t) = for *exp [i(±kr - ut)] 


V 2 / = 


1 

dt 2 


( 2 . 2 ) 


Using the Laplacian in spherical coordinates and noting that / has no dependence on 6 or 
(j) then the left hand side of Eq. 2.2 is 


7 2 r _ 1 ® f 2 


V f r 2 dr p dr ) ’ 

1 d 


= /o^^h 2 (~ r 2±r lik ) exp[i(±/cr - ut)}} , 

= /o^TT {( _1 ± rik ^> ex P[*(±^ r - wi )]} > 
r z or 

= fo ~2 {{rik ±ik — rk 2 ) exp [i(±kr — wf)]} , 


= ~k 2 f. 


(2.3) 

(2.4) 

(2.5) 

(2.6) 
(2.7) 


Now since v p = ui/k the right hand side of Eq. 2.2 is 


v 2 dt 2 


—yj = -k 2 f 


and so the 3D wave equation is satisfied. 


(2.8) 


2-2 A monochromatic plane wave E(r, t) = Eo exp[*(k • r — uit)] is travelling in the +z direction 
through a lossless linear medium with relative permittivity e r — 4 and relative permeability 
Hr « 1 and is polarised in the x direction. The frequency is v = 1 GHz and E has a 
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maximum value of +10 -3 V m 1 at t = 5 ns and z = 1 m. 

(a) Find the angular frequency, phase velocity, wavenumber, wave vector, and wavelength. 

(b) Obtain the instantaneous expression for E(r,£) valid for any position and time. 

(c) Obtain the instantaneous expression for H(r, £) valid for any position and time. 

(d) Find the Poynting vector and its time-averaged value. 

(e) Find the locations where E x is maximum when t — Os. 

Solution 

(a) The angular frequency is uo — 2nis — 2n x 10 9 = 6.24 x 10 9 rad/s. 

The phase velocity is 


1 1 c 

y/JE yj [l r \l{)E r E^ yjfl r £ r 


. 1.5 X 10 8 m/s. 


(2.9) 


The wavenumber is 


k = 


uj 

Vp 


6.238 x 10 9 
1.5 x 10 8 


41.9 rad/m. 


( 2 . 10 ) 


The wave vector has magnitude k and is in the direction of wave propagation, i.e. k = 
41.89z rad/m. 

The wavelength is A = 2n/k = v p /v = 1.5 x 10 8 /10 9 = 0.15 m. 

(b) Taking the real part of the monochromatic plane wave, and noting that the wave is 
travelling in the +z direction, has maximum Eq = 10 -3 V m -1 and is polarised in the x 
direction, 


E(r, t) = 10 -3 cos [41.89* - 6.238 x 10 9 t + 5] x V m _1 . 


( 2 . 11 ) 


To find the phase constant 5 we use the information that E is maximum when t — 5 ns 
and * = 1 m, and this occurs when the phase is 2nn (n is any integer), so 

41.89 x 1 — 6.238 x 10 9 x 5 x 10 -9 + 8 = 2nn (2.12) 


giving 5 = (2nn — 10.70) rad. The fields are unaffected by our choice of integer n. So, for 
convenience, setting n = 0 we find 

E(r, t) = 10“ 3 cos [41.89^-6.238 x 10 9 t-10.70] X Y/m. (2.13) 
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(c) We use B = k x E/ 00 — k x E/up, then 

B(r, t) = z x [lCT 3 cos (41.892 - 6.238 x 10 9 f - 10.70) x] /1.5 x 10 8 T, (2.14) 

= 6.67 x 1CT 12 cos (41.892 - 6.238 x 10 9 t - 10.70) y T. (2.15) 

Now we can readily obtain H = B /n = B///.o, 

H(r, t) = 5.31 x 10 -6 cos (41.892 - 6.238 x 10 9 i - 10.70) y A/m. (2.16) 

(d) The Poynting vector is 

S(r,i) = E x H, (2.17) 

= (10“ 3 X) x (5.31 x 10 _6 y) cos 2 (41.892 - 6.238 x 10 9 t - 10.70) , (2.18) 

= 5.31 x 10 -9 cos 2 (41.892 — 6.238 x 10 9 t — 10.70) z W/m 2 . (2.19) 
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Because the time-averaged value of cos 2 (cot) is 1/2, we have 
(S) 2.65 x 10 -9 z W/m 2 . 


( 2 . 20 ) 


(e) We need to find where E x is maximum at t = 0. This occurs where the phase is 2nn 
with n being an integer, i.e. E x is maximum where 


2nir = 41.89 xz- (6.238 x 10 9 ) x 0 - 10.70. 

.-. ^ = (10.70 + 2rwr)/41.89 = 0.255 + 0.15 n m. 

Note that the term 0.15 n is simply nX. 


( 2 . 21 ) 

( 2 . 22 ) 


2-3 Given the electric fields for the following polarisations, 

E(r,£) = (xi?o,£ + y^o , y ) cos (kz — out + 5) (linear), (2.23) 

E(r, t) — Eo[5icos(kz — cot + S x ) — ysin (kz — ut + 5^)] (left circular), (2.24) 

E(r, t) = x£^o,ic cos (kz — uot + 6 X ) + y^o,y sin (kz — ut + 6 X ) (right elliptical) (2.25) 

with Eq^ > E'o^, find the instantaneous and time-averaged energy densities and Poynting 
vectors. In each case, assume the wave is propagating in a medium with permittivity e and 
permeability fi. 

Solution 

The energy density, magnetic field, magnetic intensity and Poynting vector of an EM wave 


(2.26) 

—kxE(r,t) = v / ^ukxE(r,t), (2.27) 

Vp 

E(r, t) = ZkxE(r,i), (2.28) 

= E(r, t) x H(r, t) = [E(r, t)} x \zk x E(r, t)l = Z|E(r,t)| 2 k, (2.29) 

where v p — 1 / yfjle is the phase velocity and Z = \Jej //, is the wave impedance of the 
medium. 


are given by 


u(r ,t) 

= e|E(r,t)| 2 , 

B(r ,t) 

= ui~ l k x E(r, t) - 

H(r ,t) 

B(r, t) _ IY j 

\i V M 

S(r ,t) 

= E(r, t) x H(r, t) 


23 


Download free eBooks at bookboon.com 







Essential Electrodynamics: Solutions 


Electromagnetic waves in empty space and linear dielectrics 


(a) Linear polarisation: 


E(r, t ) = (5tE 0tX + yE 0tV ) cos (kz -ut +8), 

(2.30) 

u(r,t ) = s(Eq x + El y ) cos 2 (kz — ut + 8), 

(2.31) 

( u ) = + E l,y)i 

(2.32) 

S(r ,t) = Z(Eq x + cos 2 (kz — ut + 5) z, 

(2.33) 

(5) = l -Z(El x + El y ) z. 

(2.34) 


(b) Left circular polarisation: 

E(r, t) — Eo[5£cos(kz — uot + S x ) — y sin (kz — l ut + 5 X )\, 

u(r,t) — sEq [cos 2 (kz — ut + 5 X ) + sm 2 (kz — uot + 5 X )] — sE$, 

(u) = sEq, 

S(r ,t) = ZEqZ, 

(S) = ZE 2 0 z. 


(b) Right elliptical polarisation (E^ y > Eq^ x ): 

E(r, t) = cos (kz — uot + S x ) + yEo, y sin (kz — uot + 5 X ), 

u( r, t) = z[E^ x cos 2 (kz — uot + S x ) + E^ y sin 2 (kz — ut + 5 X )], 
u(r, t) = £ [El x + (El y - El x ) sir?(kz -ut + <5 X )] , 

( u ) = 2 £ ( E l x + E ly)i 

S(r, t) = Z [El x + (El y - El x ) sin 2 (kz - ut + $*)] % 

(S) = l -Z(El x + El y ) z. 


2-4 For a monochromatic EM plane wave incident on a plane interface between two dielectrics 
describe what is meant by perpendicular (a, s) and parallel (7r, p ) polarisation states. 
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Include a suitable diagram in your answer. 

Solution 

In transmission/reflection of light at an interface between two dielectrics, the incident ray, 
transmitted ray and reflected ray are in the plane of incidence, which is defined as the plane 
containing the incident ray and the normal to the interface. If the plane in which the elec¬ 
tric field oscillates is perpendicular to the plane of incidence, the wave has perpendicular 
(cr, s) polarisation, and if the the plane in which the electric field oscillates is parallel to 
the plane of incidence, the wave has parallel (tt, p) polarisation. 


a (s) Polarization n (p) Polarization 

incident wave incident wave 
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2-5 Define what is meant by the amplitude reflection and transmission coefficients, and by the 
reflectance and transmittance. 

Solution 

The amplitude reflection and transmission coefficients relate the complex amplitudes of 
the electric fields of the reflected and transmitted waves to the complex amplitude of the 
electric field of the incident wave. For example for perpendicular polarisation, 

= r±E ^~, Ej- = t ± E±. 

Reflectance is defined by the ratio of the the energy flux reflected by unit area of the 
interface to the energy flux incident on unit area of the interface. 

Transmittance is defined by the ratio of the the energy flux transmitted through unit area 
of the interface to the energy flux incident on unit area of the interface. 


2-6 (a) Use the amplitude reflection coefficient for parallel polarisation 


, x ( nt/rii ) cos 6i-J ( n t /n*) - sin 2 0* 

r\\ (Oi) = —j| = - ~r — 

E i (nt/rii ) 2 cos Oi + yj ( n t /rii ) 2 - sin 2 


(2.35) 


to show that Brewster’s angle is given by both 


cos 6u = — y - and tan 6U = —. 

Vi + (n t /n ? ) 2 


(2.36) 


(b) Show that 

tan 6 b — — (2.37) 

rii 

can also be derived from Snell’s law and by requiring the angle between the reflected and 
transmitted rays to be 7r/2. 
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Solution 

(a) Brewster’s angle is the angle of incidence such that r\\(0 b) — 0, so that 


(nt/rii) 2 cos 6 b = \j Ot/n*) 2 - sin 2 0 b, 
(n t /rq) 4 cos 2 6 b = ( n t /rq ) 2 - (1 - cos 2 0 B ), 


cos 2 6 b 


(nt/rii) -1 = (nt/rii) - 1, 


cos 2 9 b 


{nt/rn) 2 — 1 {nt/rii) 2 + 1 = {n t /rii) 2 - 1, 


cos 2 


(' n t/ n i ) 2 + 1 


= 1 , 


COS = 


a /1 + {nt/ni) 2 


The alternative formula we get from Eq. 2.42 above, from which 


(2.38) 

(2.39) 

(2.40) 

(2.41) 

(2.42) 

(2.43) 


cos 2 6b {nt/nt) 2 = 1 — cos 2 6b, 
cos 2 6 b ( nt/nt ) 2 = sin 2 d B , 
tands = —• 

Tli 


(2.44) 

(2.45) 

(2.46) 


(b) Eq. 2.46 can be derived from Snell’s law and requiring the angle between the reflected 
and transmitted rays to be n/2. For this case 6t = (tt/2 — 6/), and so from Snell’s law we 
find 


rii sin 6^ = rit sin 0 t , 

(2.47) 

Hi sin 6 b = n t sin^ - d B ^j , 

(2.48) 

Tii sin^ = Tit COS#B? 

(2.49) 

j . nt 

tan Ub — —• 

Ui 

(2.50) 
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2-7 Give a formula with explanations for the reflectance and transmittance in terms of the 
magnitude of the Poynting vectors of the various waves, and their angles with respect to 
the normal to the interface. How is the reflectance related to the amplitude reflection 
coefficient? How is the transmittance related to the reflectance? Why? 

Solution 

In terms of the Poynting vectors of the three waves, 

R _ s r cos e r T _ s t cos e t 

Si cos Oi ’ Si cos 6i 

Note that the area of the interface as seen by each of the three waves is the projected area, 
and the hence the relevant cos 9 factor is needed. 

The reflectance is related to the amplitude reflection coefficient by R = |r| 2 . 

The transmittance is related to the reflectance by R + T — 1 as required by energy conser¬ 
vation. 
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2-8 Explain the physical meaning of the critical angle, and derive its formula. 

Solution 

From Snell’s law rit sin 6 t = rii sin Oi 

. a n i . a 

sin u t — — sin . 
n t 

For internal reflection, Snell’s law can not be satisfied for angles of incidence greater than 
the critical angle, 

a . (nt 

u c — arcsm — 

\ n i 

for which sin 6^ > 1. For larger angles of incidence, there will be no transmitted wave, 
i.e., there will be total internal reflection. 
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3 Electromagnetic waves in dispersive media 


3-1 The resistivity of silver is p = 1.6 x 10 8 12 m, and its permeability is /j, — 0.9998/io, find 
the reflectance for light of wavelength 500 nm at normal incidence from air. 

Solution 

For perpendicular and parallel polarisations the amplitude reflection coefficients for normal 
incidence are 


r ± (0°) 


1 - nt/rij 
1 + n t /rii' 


ni(o°) 


nt/rij - 1 

rtt/rii + 1 ’ 


(3.1) 


and we see that |rj_(0°)| = |rii(0°)|. Hence, we can use either amplitude reflection coeffi¬ 
cient in calculating the reflectance, as we must be able to do since parallel and perpendicular 
polarisations are undefined for 6i = 0°. Then, for external reflection from air we have 


1 — n(uj) 1 — Re{n} — ilm{n} 
1 + n(uj) 1 + Re{n} + i Im{n} 


(3.2) 


Using the rule that the complex conjugate of a quotient of two complex numbers is the 
quotient of the complex conjugates of the two complex numbers, the reflectance is 


R( 0°) = |r(0°)| 2 


(1 — Re{n}) 2 + (Im{n}) 2 
(1 + Re{n}) 2 + (Im{n}) 2 


(3.3) 


From the wave number of a good conductor 
k(u) « (1 + i)2~^ 2 (paw) 1 / 2 


(3.4) 


we are able to estimate the refractive index, 


n(cj) ~ 



(3.5) 


A wavelength of 500 nm corresponds to an angular frequency of uo — 2n x 3 x 10 8 /500 x 10 9 
3.77 x 10 15 rad s -1 . We see that for silver 


Re{n} « I m{n} = 3 x 10 8 


1 0.9998 x (4tt x lO" 7 ) x (1/1.6 x 10" 8 ) 


2 x (3.77 x 10 15 ) 


= 30.6. 


(3.6) 
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Hence, using Eq. 3.3 the reflectance is 


R( 0°) 


(1 — Re{n}) 2 + (I m{n}) 2 
(1 + Re{n}) 2 + (Im{n}) 2 


(3.7) 


3-2 Show that the time-averaged power density (E • J) in a dilute plasma is zero. 

Solution 

The conductivity of a dilute plasma is purely imaginary, 

n e e 2 

<7 = l , 
m e (jj 

so that the current density and electric field are 90° out of phase. Consequently, the 
time-averaged power will be zero, as shown below 



E J = RejEoe-^ 4 } • Re {e^ /2 |cr| E 0 e“ iwt } , 

(3.9) 

= {E 0 cos ■ | |cr| E 0 cos (| - tot) |, 

(3.10) 

= \ct\Eq cos(u;t) sin (cut) , 

(3.11) 

E • J = - \<j Eq sin(2c ut) 

(3.12) 

;e-j) = o. 

(3.13) 


Hence, there is no resistive energy loss. 


3-3 At midnight and noon, at a certain location and date, the electron number density in the 
ionosphere was as given in the plot below. 

(a) Label the x-axis at the top of the plot in terms of the plasma frequency v v (MHz). 

(b) Find the height of the reflecting layer at 1 MHz and 3 MHz at midnight and noon. 

(c) Find the minimum frequencies for communication with orbiting satellites at midnight 
and noon. 

(d) Estimate the total electron content in (electrons m -2 ) at midnight and noon. 

(e) An instantaneous broad-band pulse of radio-frequency interference (RFI) is emitted 
overhead by an orbiting satellite and observed by an terrestrial detector with a bandwidth 
from 1.2—1.8 GHz. What is the duration of the observed pulse? 
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log[n e /m' 3 ] 


Solution 

(a) The plasma frequency in Hz is v v — oj p /2n 

_ 1 f n e e 2 \ 1/2 
Up ~ 2tt \eomJ 


(3-14) 
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and so for n e = 10 9 m -3 the plasma frequency is 2.8 x 10 5 Hz and log(z/ p /MHz) = —0.55. 
Because v p oc n~J 2 two decades of n e correspond to one decade of i/ p , so we now have the 
information needed to relabel the plot as shown below. 


I°g( v p / MHz) 

-0.5 0.0 0.5 1.0 



9 10 11 12 13 

log[/i e / m" 3 ] 


(b) Reading off the plot, we obtain the height of the reflecting layer shown in the table 
below. 


Time 

1 MHz 

3 MHz 

midnight 

noon 

240 km 

80 km 

300 km 

85 km 


(c) Reading off the plasma frequencies at the maximum n e values, the minimum frequencies 
to avoid reflection by the ionosphere are 6.3 MHz (midnight) and 18 MHz (noon). 

(d) The total electron content (TEC) is f n e (h)dh , which we can estimate by tabulating 
n e (in units of 10 12 m -3 ) at every 100 km, and multiplying the sum by Ah = 10 5 m as 
follows. 
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h (km) 

noon 

midnight 

100 

1.4 


200 

3.9 


300 

3.6 

0.17 

400 

1.8 

0.50 

500 

0.8 

0.34 

600 

0.4 

0.22 

700 

0.3 

0.16 

800 

0.3 

0.14 

900 

0.2 

0.13 

Sum 

12.7 

1.66 


Hence, the TEC at noon is (12.7 x 10 12 ) x 10 5 = 1.27 x 10 18 electrons m 2 , whereas at 
midnight the TEC is 1.66 x 10 17 electrons m -2 . 

(e) Using the result for a given dispersion measure, which is defined in a similar way to the 
TEC, we have the arrival time of the pulse at a given frequency 


p 2 

t(v) = —a -x TEC x v~ 2 = (1.34 x 10“ 7 ) x TEC i/“ 2 (3.15) 

87T 2 £om e C 

such that the spread in arrival times is 

At = (1.34 x 10“ 7 )(i/f 2 - z/ 2 -2 ) x TEC (3.16) 

= (1.34 x 10- 7 ) x [(1.2 x 10 9 )- 2 - (1.8 x 10 9 )" 2 ] x TEC (3.17) 

= 2.46 x 10- 26 x TEC. (3.18) 


Hence, the pulse is spread out in time by 3.13 x 10 8 s (noon) or 4.09 x 10 9 s (midnight). 


3-4 For pulsar PSR J1817-3837 the observed position angle (PA) of linear polarisation with 
respect to the North Celestial Pole is plotted below vs. frequency. The dispersion measure 
is 102.85 pc cm -3 , where 1 pc=3.09 x 10 16 m (1 parsec) is the distance unit used by 
astronomers. Find the average value of the parallel component of the magnetic field along 
the line of sight to this pulsar. 
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PSRJ1817-3837 



Solution 

The first step is to determine the rotation measure for this pulsar using the data given in 
the plot. Even though we expect a v~ * 1 2 dependence of position angle, the range of frequency 
is sufficiently small that we can make a linear fit (shown below). 
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1250 1300 1350 1400 1450 1500 

v (MHz) 

Reading off data at two frequencies we can tabulate it together with wavelength as shown. 


V (MHz) 

A (m) 

A 2 (m 2 ) 

PA (degrees) 

1500 

0.2 

0.04 

-40 

1250 

0.24 

0.0576 

+64 


The observed rotation measure is 
A PA 104° 

RM — - ^ 5909° m 2 103radm 2 . (3.19) 

A A 2 0.0176 m 2 v ; 

The rotation measure is related to the magnetic field and electron density along the line of 
sight to the pulsar by 


RM 


2.63 x 10 -13 


J n e (r) B(r) • dr 


(rad m 2 ). 


(3.20) 


Hence, for PSR J1817-3837 we must have 


J n c (r)B(r) • dr = RM/( 2.63 x 10“ 13 ) = 3.92 x 10 14 T m” 2 . 


The observed dispersion measure is 


(3.21) 


DM 


J n e (r) dr = 102.85 x (3.09 x 10 16 ) x 10 6 


3.18 x 10 24 m -2 . 


(3.22) 
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Hence, the average value of the parallel component of the magnetic field along the line of 
sight to this pulsar is 


B J n jr)B(r)-dr 

J n e (r) dr 


1.23 x 10" 10 T. 


(3.23) 


3-5 We are receiving radio signals from an interplanetary space probe which is far from the 
Sun. It is currently on a trajectory such that it has been eclipsed by the Sun, but it is now 
emerging from behind the Sun. The space probe regularly broadcasts instantaneous broad¬ 
band pulses of radio-frequency emission every second, and it also broadcasts spacecraft 
experiment instrument data over a narrow frequency band of width 100 kHz extending from 
10,000.0 MHz to 10,000.1 MHz. It transmits using a dipole antenna aligned perpendicular 
to the plane of the ecliptic (the plane containing the Earth’s orbit). We are just renewing 
radio contact as the space probe is emerging from behind the solar limb. 

The Sun’s corona is an extremely hot (> 10 6 K) plasma which is highly variable, and 
has dynamic coronal loops of magnetised plasma. The corona is located above the Sun’s 
photosphere. At the time we are observing, the base of the corona in front of the space probe 
as it emerges from behind the Sun has an electron number density of n e (i?© ) - 3xl0 15 m -3 . 
The corona’s a scale height is about H ~ 10 8 m, such that the electron number density 
decreases with height as n e {r ) = n e (i?©) exp[—(r — Rq)/H]. Imagine that, at the time we 
are observing, the Sun’s magnetic dipole moment is m = 2 x 10 29 A m 2 and it happens to 
be pointing directly towards Earth. 

(a) What is the minimum radio frequency that we are able to receive from the space probe? 

(b) What is the additional delay of the pulses due to propagation through the corona? 

(c) By how much has the plane of polarisation of the EM wave carrying the data signal 
rotated due to propagation through the corona? 

[Make what you think are reasonable approximations - do not attempt a rigorous calcula¬ 
tion. The radius of the solar photosphere is i? 0 = 6.955 x 10 8 m.] 

Solution 

The geometry is as shown below. 
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(a) The plasma frequency is 

(3. 

and so at the base of the corona where n e (i?©) = 3 x 10 15 m -3 the plasma frequency is 
3.1 GHz. No communication can take place below this frequency. 

(b) To estimate the dispersion measure we could use the maximum electron number density 
and a path-length of 2 L corresponding to the path through the corona below height H (see 
diagram above), where 
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L = ^J(Rq + H f - R% = 3.9 x 10 8 m. (3.25) 

Hence, the dispersion measure is DM ~ 2Ln e (i?©) ~ 2.3 x 10 24 m -2 . 

Using the standard result for the time delay for a given dispersion measure, we have the 
arrival time of the pulse at different frequencies 

e 2 

£( v) = 8yr 2 g m c x DM x = (1.34 x 1CT 7 ) x DMU 2 (3.26) 

Over the bandwidth of the detector the delays are 

£(10,000.0 MHz) = (1.34 x 10 -7 ) x (2.3 x 10 24 ) x (1.0000 x 10 10 ) -2 = 3.0820xl(T 3 s, 
£(10,000.1 MHz) = (1.34 x 1(T 7 ) x (2.3 x 10 24 ) x (1.0001 x 10 10 ) -2 = 3.0814x 1CT 3 s. 


Hence, the start of the observed pulse is delayed by ^ 3 ms and the duration of the pulse 
is ^ 600 /is. 

(c) To find the angle of rotation of the plane of polarisation, we need to estimate the parallel 
component of magnetic field by, say, the surface magnetic field at the magnetic equator 
r = (i? 0 , 90°, </>), where the pole of these spherical coordinates corresponds to the magnetic 
axis of the Sun. Then, we can obtain the surface field from the magnetic dipole moment 
m Q using 


B(r) 


/xo 3r(m 0 • r) - r 2 m Q 
47T r 5 


(3.27) 


and so at the Sun’s magnetic equator the surface magnetic field magnitude is 


B = 


Mo rap 
4t r R% 


4tt x 10“ 7 2 x 10 29 

47t (6.996 x 10 8 ) 3 


= 5.8 x 10" 5 T. 


(3.28) 


The rotation measure is related to the magnetic field and electron density along the line of 
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sight to the pulsar by 


RM 


2.63 x 10" 13 

j n e ( r) B(r) • dr (rad m 2 ) 

(3.29) 

(2.63 x 10“ 13 ) 

* x (3 x 10 15 ) x (5.8 x 10- 5 ) x (2 x 3.9 x 10 8 ) 

(3.30) 

3.6 x 10 7 rad 

m 2 . 

(3.31) 


The rotation measure is defined through 
A PA 

RM = (3.32) 

Hence, the plane of polarisation will have rotated by 

/q x 10 8 i \ 2 

PA(10,000.0 MHz) = (3.6 x 10 7 ) x f ’ J = 32351.522 rad, (3.33) 

/ 3 x 10 8 i \ 2 

PA(10,000.1 MHz) = (3.6 x 10 7 ) x ——-A— = 32345.053 rad. (3.34) 

\ 1.UUU1 X 1U / 


So the plane of polarisation will have rotated by an enormous angle of the order of 30,000 
radians, and we see that over the 100 kHz receiver bandwidth the position angle changes 
by 5 radians. These numbers are of course only order of magnitude estimates, but 
they show that trying to align the a dipole antenna with the plane of polarisation makes 
no sense, and that to obtain all the signal it would be necessary to have two orthogonal 
polarisations and appropriately combine the two signals. 
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4 Waveguides 


4-1 A standard waveguide for the “X-band” (8.2-12.4 GHz) has internal cross section 2.286 cm 
by 1.143 cm. Find the the first two cut-off frequencies, as these will give the frequency 
range for which only the TEio mode will propagate. 

Solution 

Since a — 2b, these will correspond to c^io — cn/a and c^oi = ^20 — 2c7r/a = 2cjiq. To get 
the frequency in Hz, 

Z'lo = c/ (2a) (4.1) 

= (3 x 10 8 )/(2 x 2.286 x 10 -2 ) (4.2) 

v 10 = 6.56 GHz, (4.3) 

.-. v Q1 = 13.1 GHz. (4.4) 


4-2 Find the electromagnetic field energy per unit length of waveguide for TEio mode. Express 
your result in terms of ateo^max- Start with the fields for the TEio mode 


( 7VX \ 

—) sin (kz - ut), 


_ k . /7TX\ 

B x = -E max -srn —) sin (kz-ut), 
cj V a J 

7T / 7 TX \ 

B z = — E max —cos ( — ) cos (kz — ut). 
aid \ a J 


where k = ^/cd 2 /c 2 — ( 7 r/a) 2 . 


(4.5) 

(4.6) 

(4.7) 


41 


Download free eBooks at bookboon.com 





Essential Electrodynamics: Solutions 


Waveguides 


Solution 

The time-averaged energy in the EM field per unit length of the waveguide is 



/ J 0 {([- B ^( r > t )] 2 ) + {[B z (r,t)} 2 )} dydx, 


(4.8) 


n b 


^ sin 2 ^ sin 2 {kz — Lot)^ dy dx 

sin 2 (fez — cat)^ 



(4.9) 


Now, we can use the following: 


(sin 2 (fc 2 : — ut)) = (cos 2 (kz — ut)) — [ sin 2 f—i dx = f cos 2 f—) dx = —. (4.10) 

2 Jo \ a / Jo \ a / 2 


Then 


U e 0 ab 2 1 ab k 2 2 1 ab it 2 2 

L 2 4 max + 2 M0 4 W 2 “ aX + ^ 4 ^^max, 


(4.11) 


and using the dispersion relation given we can substitute for A: above to get 


U_ 

1 


£^£2 , 
2 ^ max 1 


1 ab f 1 7r 2 ^ ^2 , 1 ab n 2 

2iTo^W~^ 2 ) max + 2^T^2 max> 


gp ab 2 
2 4 max ' 2/x 0 4 c 2 ~ max ' 


1 ab 1 o 

+ z- 


(4.12) 

(4.13) 


Hence, remembering that /Togo = 1/c 2 , the energy per unit length of waveguide for TEio 
mode is 


U = £oab 2 
L ~ 4 £ ' max ' 


(4.14) 
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4-3 Consider a hollow standard WR-159 F-band waveguide which has an internal cross section 
of 40.386 mm x 20.193 mm. 

(a) A 2 GHz signal is fed into the waveguide. Will this cause a wave to propagate along 
the waveguide? 

(b) Assuming the wave would be in the TEio mode, what is the wavenumber? Discuss the 
implications of your answer. 

Solution 

(a) The lowest cutoff frequency is for the TEio with cjiq = cn/a = 2.33 x 10 10 Hz, hence 
v io = c/2a = 3.71 GHz. 2 GHz is below the cut-off frequency for the TEio mode, and so 
the wave will decay exponentially with distance. 

(b) The wavenumber is obtained from the dispersion relation with uo = 2n x 2 rad s -1 

k — -yj a; 2 — — ±i65.6 m _1 . (4.15) 

Because this is purely imaginary the wave will not propagate, and the electromagnetic field 
will decay exponentially as e ~ lrn { k } z if the wave was launched in the +z direction. The 
field attenuation length is 1/1 mn{k} — 15.3 mm. 
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4-4 Consider the case of launching an EM wave with frequency uo < cjio into a rectangular 
waveguide. Assuming the wave will attempt to propagate as a TEio wave travelling in 
along the waveguide (in the z direction), 

(a) Derive all the components of the electromagnetic field for the case where k is purely 
imaginary. 

(b) Obtain the Poynting vector and discuss the energy flow. 

Solution 

(a) We will need to obtain the (real) EM field E(r, t ) and B(r, t ) of the TEio mode for the 
case where k is purely imaginary. We know that for all the TE modes E®(x,y) = 0 and 
that the amplitude function of the longitudinal component of the magnetic field for the 
TEio mode is 



(4.16) 


We can get all the other field components from B® and E® using 



(4.17) 



(4.18) 



(4.19) 



(4.20) 


where we have used 


^ (jo 2 — k 2 c 2 





(4.21) 


We need to multiply these field amplitude functions by exp [i (kioz(u) — cut)] to get the 
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(complex) EM field of a TEio wave in the ^-direction: 


E x (r,t) = 0, (4.22) 

E y (r,t) = + Bioioj 0 sin (0 e ^ kz ~ ut \ (4.23) 

E z (r,t) = 0, (4.24) 

B x (r,t) = -B 10 ik 0 sin (0 e i{ - kz ~^\ (4.25) 

B y ( r ,t) = 0, (4.26) 

B z ( r,i) = + Bio cos (0 e’( fc2 -^), (4.27) 

with k = y/u 2 /c 2 — ( 7 r/a) 2 . (4.28) 


In the present problem, fc(a;) = zlm{fc(a;)}, and so the (real) EM field of the TEio wave in 


the z-direction is 

£k(r,i) = 0, (4.29) 

E y (r, t) = Re |+ £>io iu ^ sin ^ , (4.30) 

= - Blow 0 sin (0 Sin(-wt), (4.31) 

= B max sin (~) e“ Imffc}2 sin(wt), (4.32) 

E z (r,t) = 0, (4.33) 

B x (r,t) = Re^-B w ik (^) sin (^) e* (fe2_wt) } , (4.34) 

= +Bi 0 Im{fc} 0 sin (0 e“ Im W 2 sin (-cot), (4.35) 

= —E max — sin (—) e -Mfe}« sin(wt), (4.36) 

(jJ V CL / 

B y (r,t) = 0, (4.37) 

B*(r,*) = Re {+Bi 0 cos (0 e^ 2 ”^} , (4.38) 

= +Bio — cos (—) e _Im i fc } 2 cos(-wt), (4.39) 

7r V a J 

= -^max— cos (—) cos(wf). (4.40) 

cja V a / 
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(b) The Poynting vector is S = E x H = ( E y H z St — E y H x z) for this problem, then 


E y (r,t) = £ max sin(^) e lm ^ z sin(ut), (4.41) 

H x (r,t) = — E'max —1-1 sin e~ lm ^ z sin (ut), (4.42) 

cj/io V a J 

H z (r,t) = -E m ax - cos (—) e -M^}2 C os (cot). (4.43) 

bj/iQa V a / 

Hence, the non-zero components of the Poynting vector are 

<Sz(r, t) = ^Lx^^ sin 2 (—) e - 2I -W-i sin(2o;i). (4.44) 

S*(r,t) = e - 2 ta{^ 1 in(2ajt)> (4.45) 

(jj fjjQd A \ a J A 

where we have used cos H sin A = sin(2A)/2. 


The result is that S oscillates at frequency 2cj and decays with increasing z, and (S) = 0. 
The intensity attenuation length is (21 m{k})~ 1 . 
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5 Radiation and scattering 


5-1 A walkie-talkie has a short centre-fed dipole antenna 12 cm long and radiates 3 W at 
450 MHz. 

(a) What is the radiation resistance of the antenna? Assume the antenna has been balanced 
by a suitable inductor such that its impedance is purely resistive, and that the transmission 
line from the transceiver is matched to the impedance of the antenna. 

(b) What is the value of the peak current in the antenna? 

(c) What is the peak electric field at the location of a receiving antenna at distance r? 

(d) With the transmitting walkie-talkie held vertically, an identical receiving walkie-talkie 
also held vertically at a distance of 2 km away can barely receive the transmission, i.e. 
the range is 2 km. What would be the range if: (i) the transmitting antenna were tilted 
at 60° towards the receiving antenna, (ii) the transmitting antenna was vertical but the 
receiving antenna was tilted at 60° towards the transmitting antenna, (iii) the transmitting 
antenna was vertical but the receiving antenna was tilted at 60° to the vertical in a plane 
perpendicular to the line of sight to the transmitting walkie-talkie? 

Solution 

(a) The wavelength is A = 3 x 10 8 /4.5 x 10 8 = 0.667 m, giving d/X = 0.18 and so we 
shall use formulae relevant to a short (d <C A) centre-fed dipole antenna. The radiation 
resistance of this short centre-fed dipole antenna is then approximately 


^rad = 197 



= 197 x 0.18 2 


6.3811. 


(5.1) 


(b) We can use the time-averaged power (. P) = 3 W, together with the radiation resistance, 
to obtain the peak current Iq as follows 


II = 2 (I 2 ) = 2 


( P) 

-Rrad’ 


h = 


2x3 

6.38 


= 0.97 A. 


(5.2) 


(c) The antenna will emit electric dipole radiation, and for this short centre-fed dipole 
antenna the dipole’s amplitude is 


Po 


Iod 
2 u 


0.97 x 0.12 

2 x (2 tt x 4.5 x 10 8 ) 


= 2.06 x 10“ n C m. 


(5.3) 
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For electric dipole radiation the electric field, and field amplitude, are 


E(r ,t) = - 


^ 2 \ e i(kr—ut) 


47T£o ) 


PO 


sin ee. 


En = 


( u- \ Po . 

I --ft I — sin0 0, 


\47T£o c 2 J r 

(27r x 4.5 x 10 8 ) 2 

_ 4tt x (8.85 x 10- 12 ) x (3 x 10 8 ) 2 

-16.5^9 (VO. 
r 


2.06 xlO- 11 ^?. 


(5.4) 

(5.5) 

(5.6) 

(5.7) 


(d) When both antennas are vertical, in a horizontal plane, and 2 km apart, the electric 
field component parallel to the receiving antenna is the minimum detectable, 

ii 16 5 

E l L = sin 90° = 8.25 x 10“ 3 V m" 1 . (5.8) 

(i) For the transmitting antenna tilted by 60° towards the receiving antenna we have 
9 = 60° in Eq. 5.7, and 6 is pointing vertically down (parallel to the receiving dipole 
antenna). Hence, at the detection threshold we must have 


ii sin 60° 

^iin = 16 - 5 ^—, ••• r = 1.73 km. (5.9) 

(ii) If the transmitting antenna is vertical and the receiving antenna is tilted by 60° towards 
the transmitting antenna, the electric field is vertical but the component parallel to the 
antenna is reduced by cos 60°. Hence 


<in = 16-5 


sin 90° cos 60° 


r = 1.00 km. 


(5.10) 


(iii) If the transmitting antenna is vertical and the receiving antenna is tilted by 60° from 
the vertical perpendicular the direction to the transmitting antenna, the electric field is 
vertical but the component parallel to the antenna is reduced to F?ocos60°. Hence, as 
before, 


Kin = 16-5 


sin 90° cos 60° 


r = 1.00 km. 


(5.11) 
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5-2 Starting with the Clausius-Mossotti formula and the time-averaged dipole radiation power 

“■» = <p > = (5 - i2) 

find the cross section for scattering of light at A 500 x 500 nm wavelength by air molecules, 
and its associated mean free path l/(n v a). Dry air at STP has p = 1.30 kg m -3 , and 
n — 1.00029, and its molecular weight is 29. 

Solution 

A molecule acquires an electric dipole moment in the presence of an electric field p = <a m E. 
For the case of a dilute gas the Clausius-Mossotti equation yields 

= — (e r - 1) = -(n 2 - 1) (5.13) 

Ti v Ti v 

where n v is the number density of molecules, e r the relative permittivity and n is the 
refractive index. 



power 


By 2020, wind could provide one-tenth of our planet’s 
electricity needs. Already today, SKF’s innovative know¬ 
how is crucial to running a large proportion of the 
world’s wind turbines. 

Up to 25 1 of the generating costs relate to mainte¬ 
nance. These can be reduced dramatically thanks to our 
^sterns for on-line condition monitoring and automatic 
lul^kation. We help make it more economical to create 
cleanSbjheaper energy out of thin air. 

By snSfinoi our experience, expertise, and creativity, 
industries can boost performance beyond expectations. 
1 Therefore we need the best employees who can 
fcneet t.hjs challenge! 


The Power of Knowledge Engineering 


Brain 


Plug into The Power of Knowledge Engineering. 
Visit us at www.skf.com/knowledge 




Click on the ad to read more 


Download free eBooks at bookboon.com 




Essential Electrodynamics: Solutions 


Radiation and scattering 


In the presence of an oscillating electric field due to an incident monochromatic plane 
wave E 0 cos(k • r — uot) the molecule becomes an oscillating dipole which radiates with 
time-averaged power 


(■ p) 


127T£()C 3 


/ g 0 \ 2 (n 2 - l) 2 Ega> 4 
\n v J 127T£oc 3 


(5.14) 


Dividing ( P) by the time-averaged energy flux of the incident wave ( S) — we 

obtain the cross section 


^mol 


87T 3 (n 2 - l) 2 
3 n 2 


A” 4 . 


(5.15) 


From the data on dry air given, 


(n 2 - l) 2 
n v 

cr(A) 


= 3.4 x 10" 7 , 


p _ 1.30 

29 u ~ 29 x 1.66 x 10 “ 27 

. 8 ^( n 2 - l )\_ 4 


= 2.6 x 10 25 m -3 . 


-A - 


3 n 2 v 
8tt 3 3.4 x 10 


3 (2.6 x 10 25 ) 2 


-7 


x (500 x 10- 9 A 500 )- 4 , 


= 6.65 x 10 31 A 50 q m 2 . 


(5.16) 

(5.17) 

(5.18) 

(5.19) 

(5.20) 


The mean free path is l/n v a = 58A5 00 km, where A 500 = A/(500 nm). 


5-3 In the presence of an applied uniform electric field E = Eqz the radial component of the 
electric field outside a perfectly-conducting sphere of radius a centred at the origin becomes 


E r (r,6,<t>) = E 0 



COS#. 


(5.21) 


Find the cross section for scattering of monochromatic EM waves by a perfectly-conducting 
sphere of radius: (a) a « A, and (b) a^> A. Assume electric dipole radiation is responsible. 


Solution 

(a) For a monochromatic plane wave with electric field E(r,£) = Eqcos(A:x — cjt)z with 
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k 1 a, the applied electric field over the entire sphere can be approximated by E(a, 0, </>, t) ~ 
cos(cc;t)z. 

Then, using Gauss’ law, the surface charge density is 


<r(a, 0, 0, £) = £oE(a, 0, 0, £) • ? = £oE^a, 0,0,£) = 3^o^o cos(c<;£) cos 0. (5.22) 

The dipole moment will be in the z direction, and its z component will be 


Pz(t) 

= z • j) r <r(r, t) dS , 

(5.23) 


= cos (cot) J (a cos 9)(3eoEo cos 9) 2na 2 d(cos0), 

(5.24) 

Pz(t ) 

= 67ra 3 £oEo cos(ujt) J cos 2 9 d (cos 0). 

(5.25) 

P z(t) 

= 47ra 3 6:oEo cos(c ut) z. 

(5.26) 


Using Larmor’s formula for electric dipole radiation, the time-averaged power is 


(p) 


12'keqc 3 


(47ra 3 £o) 2 El co 4 Aeo^eqEI cu 4 

127T£oc 3 3c 3 


(5.27) 


Because the dipole radiation pattern is proportional to sin 2 0 where 0 is the angle to the 
acceleration direction, the radiation will be zero in the direction of acceleration, i.e. in the 
zb z direction and maximum at all directions in the x—y plane. 

The time-averaged Poynting vector for the EM wave is 


(S( r )) = ^£oEqCZ. 


(5.28) 


Hence, the scattering cross section due to electric dipole radiation is 


EDR 




(_p± 

(S) 


/ 47ra 6 eo£’o 

[ 3^ 


-£oE 0 c ) = 


= ^(*V)a 2 . 


(5.29) 


Since the wavenumber has units (m _1 ) and a has units (m), the cross section has units 
(m 2 ) as required. The scattering cross section is proportional to u 4 which is the same 
dependence as for scattering by molecules which was found by Rayleigh. 
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Note that the cross section we have calculated is actually incomplete as it neglects magnetic 
dipole radiation from oscillating surface current loops induced by the changing magnetic 
field of the incident electromagnetic wave. Including magnetic dipole radiation would 
increase the calculated cross section by 25%. 

(b) The cross section for the case where <C a, is for wavelengths A a, in which 
case the perfectly conducting sphere may be considered a perfect reflector blocking out its 
cross-sectional area a = no 2 . 


5-4 In the presence of an applied uniform magnetic field B = Bqz a perfectly-conducting sphere 
(radius a) acquires a magnetic dipole moment mo- 

(a) What are the boundary conditions at the surface of the sphere? Find the magnetic 
dipole moment that, together with the applied magnetic field, satisfies the boundary con¬ 
ditions. 

(b) Find the contribution of magnetic dipole radiation to the cross section for scattering 
of monochromatic EM waves by a perfectly-conducting sphere of radius a<C A. 

Solution 

(a) Inside a perfect conductor E = 0 and B = 0, so that the boundary conditions on the 
electromagnetic fields at the conductor’s surface just outside the conductor are simplified 
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to E|| = 0 and = 0. That is, the electric field at the surface is normal to the surface, 
and the magnetic field at the surface is parallel to the surface. 

(b) Adding the field of a magnetic dipole at the origin with moment m = moz to the 
constant field £>qz, 


B(r) = B 0 z+^- 


3 r mo cos 9 — mgz 


= £>o(cos0r — sin 9 6) + 


Mo m o 

47r 


3 r cos 9 — (cos 9 r — sin 9 6 


At the surface of the conductor B r (a , 0, </>) must be zero, and so 


O n , W m 0 o a n 

B o cos 9 + -— t 2 cos 9 — 0. 
Ana 6 


27ra 3 _ 

m 0 =- B 0 z. 

Mo 


(5.30) 

(5.31) 


(5.32) 

(5.33) 


For a monochromatic plane wave with magnetic field B(r,£) = Bocos(kx — cot)z with 
k~ l a, the applied magnetic field near the conductor can be approximated by B(r,£) ~ 
Bq cos(c ut) z. Hence, 


m (£) — mo cos(cct) z. 


(5.34) 


For an oscillating magnetic dipole, the time-averaged radiation power is 

, , Mow 4 mg Mow 4 47r 2 a 6 2 

127TC 3 127TC 3 /ig ° 

The time-averaged Poynting vector for the EM wave is 
<S(r)> = 2 ^ B » cS ' 


(5.35) 


(5.36) 


Hence, the contribution of magnetic dipole radiation to the scattering cross section is 


MDR 

cr ar . 


(p± 

(5) 


( Mow 4 4 t r 2 a 6 d2 ^ 
Uttc 3 B °) 



2l r 

T 


(fc 4 a 4 )a 2 . 


(5.37) 


Adding this to the contribution of electric dipole radiation to the scattering cross section 
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(Exercise 5-3) we obtain 


tot MDR | EDR 

°sc = °sc +°sc 


107T 

^r 


(k A a A )a 2 . 


(5.38) 


5-5 In bremsstrahlung an energetic electron is accelerated in the electric field of an atomic 
nucleus and emits electromagnetic radiation. Consider the case of an electron of kinetic 
energy 10 keV travelling on a trajectory parallel to the z-axis, and that there is a lead 
(Pb) nucleus located at impact parameter b — 10 -10 m (distance of closest approach to 
the initial trajectory). Assume that for this impact parameter the electron’s deflection is 
small, and that a calculation of the energy lost to electromagnetic radiation can be done 
classically. 

(a) By comparing b with the de Broglie wavelength h/mv , reassure yourself that a classical 
calculation is valid for this impact parameter and velocity. 

(b) By checking whether or not the deflection angle is small, reassure yourself that approx¬ 
imating the trajectory by a straight-line is valid. 

(c) Estimate the total energy radiated by the electron. [Neglect atomic electrons.] 

Solution 

(a) From the kinetic energy, Ek = (1/2 )mv 2 , we first find the electron’s velocity, and then 
the de Broglie wavelength 


v = 


2 E k 


m P 


2 x (10 4 x 1.6 x lO" 19 ) - , 

v L = 5.93 x 10 7 m s -1 , 


= — = 


9.11 x 10 -31 
6.62 x 10“ 34 


mv (9.11 x 10" 31 ) x (5.93 x 10 7 ) 


= 1.23 


1CT 11 m. 


(5.39) 

(5.40) 


The impact parameter b = 10 10 m is about a factor of 10 larger than the de Broglie 
wavelength, and so use of classical physics is (marginally) justified. 


vt 




E 


v 


z 


(b) The straight line approximation is only valid when the impulse is small compared to 
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the initial momentum, i.e. A p± <C m e v where 

Ze 2 

AP1 = /_ co 4 TCo (i>= + t ,V) COS[e< ‘ )1 *' 

r°° z^_ b 
J-O o 47T£o (b 2 + V 2 t 2 ) 3 / 2 

Ze 2 
27T£o bv 

Since the straight line approximation is valid when Ap^ <C m e ?;, we require 
Ze 2 

b » -tt = 1.18 x 1CT 11 m 

27rsom e v z 

which is marginally satisfied since b — 10 -10 m. 

(c) Larmor’s formula for the instantaneous radiation power is 


m 


a(t) 2 e 2 

6tT£()C 3 ’ 


(5.41) 

(5.42) 

(5.43) 


(5.44) 


(5.45) 



READ MORE & 


qaiteve 


Challenge the way we run 


EXPERIENCE THE POWER OF 
FULL ENGAGEMENT... 


RUN FASTER. 

RUN LONGER.. 
RUN EASIER... 


ORDER TODAY 


PRE 



Click on the ad to read more 


Download free eBooks at bookboon.com 














Essential Electrodynamics: Solutions 


Radiation and scattering 


and so the total energy radiated by the electron as it travels from z = —oo to z = +oo is 


W = 


6tT£()C 3 J_ 


/ oo 

a(t) 2 dt, 

-oo 


/ 


/ (Ze) x e 


6tvsoc 3 J _ 00 \4:7T£o(b 2 + v 2 t 2 )m ( 

e 2 / Ze 2 x 2 
67T£()C 3 \4:7T£om t 


L 


1 


-oo 0 2 + v 2 t 2 ) 2 


dt. 


Z 2 e 6 


967r 3 £QC 3 m 2 


bt 


Z z e 


2 „6 


2b 3 (b 2 + ?; 2 t 2 ) 


7T 


+ 


tan 1 (?;£/&) 
2b 3 v 


967r 3 6:QC 3 m 2 2b 3 v 5 


2 „6 




1927^ 2 6:QC 3 m 2 6 3, ^;’ 
= 4.04 x 10- 5 eV. 


(5.46) 

(5.47) 

(5.48) 

(5.49) 

(5.50) 

(5.51) 

(5.52) 


5-6 A pulsar is a rapidly-spinning magnetised neutron star observed at radio and sometimes 
higher (optical, X-ray and gamma-ray) frequencies. Neutron stars have masses of typically 
M ^1.4 M 0 (solar mass 1 M 0 w 2 x 10 30 kg), radius R zz 10 km, moment of inertia 
ImoI — | MR 2 « 10 38 kg m 2 and a wide range of surface magnetic fields. Outside the 
neutron star the magnetic field may be approximated by that of a magnetic dipole at its 
centre making angle a to the spin axis. Pulsars typically emit one or more narrow pulses 
at the same time into each spin period P, and measurements at different epochs generally 
show that the spin is slowing down as characterised by the time-derivative of the period P. 

(a) Derive formulae in terms of P and P for the rotational kinetic energy E r Q t, the rate 
of loss of rotational kinetic energy E rot , the pulsar’s characteristic age r = E TOt /E TOt and, 
assuming the slow-down is due to conversion of rotational kinetic energy to magnetic dipole 
radiation, the pulsar’s minimum equatorial surface magnetic field B m [ n . 

(b) The radio pulsar PSR J0157+6212 has P = 2.355 s and P = 1.89 x 10 -13 (“The 
Australia Telescope National Facility Pulsar Catalogue”, Manchester, R.N., et ah, 2005, 
Astron. J., 129, 1993, and http://www.atnf.csiro.au/research/pulsar/psrcat/). Find E r 0 t, 

Eroti ~ -^rot/'E'rot and Pm in- 
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Solution 

(a) The rotational kinetic energy is 
-E'rot 2 = 2 1 -MoI 



27 T 2 lMolP 2 - 


The rate of loss of rotational kinetic energy is 


E rot = 27t 2 /m 0 /(-2)P 3 P = -47t 2 /mo/-^3 ‘ 

The pulsar’s characteristic age is 

27 T^ImqiP^ = P_ 

4vr HmoiPP- 3 2 P ' 


(5.53) 


(5.54) 


(5.55) 


For the pulsar, we have a rotating magnetic dipole of effective moment m o sum, where a is 
the angle between the pulsar’s magnetic moment and its spin axis as it is only the moment’s 
component perpendicular to the spin axis which is relevant for radiation. Furthermore, a 
rotating magnetic dipole moment is equivalent to two orthogonal oscillating magnetic dipole 
moments out of phase by 90°. 

We need to write the magnetic moment in terms of the surface magnetic field, and for a 
magnetic dipole of moment mo 

B « = £ 


3r(mo • r) — r^mo 


(5.56) 


and so at the magnetic equator of a neutron star (mo • r = 0 in the formula above) the 
surface magnetic field has magnitude 


B = 


Mo rnp 

Ai r R 3 ’ 


rn 0 


4:7V 

Mo 


R 6 B. 


(5.57) 


The time-averaged magnetic dipole radiation power of a rotating dipole is twice the power 
of a single oscillating magnetic dipole 

(ft.DE) = 2 X ^27 <5 ' 58> 

Substituting for the magnetic dipole moment in terms of the polar surface magnetic field 
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we get 


(^mdr) = 


87T 3 Ro /47T o ^ 

pu ' — R 3 B sin a 


3 c 3 P 4 \/io 


Equating this to the rate of loss of rotational kinetic energy is 


8ttVo ( 4vr 3 


3c 3 P 4 \/iq 


—P Pshm = 4tt ImoIty^ 


B 2 sin 2 <a = 47 t 2 ImoI 


P 3' 

P 3c 3 P 4 /i 2 


B sin ct — 


P 3 87T 3 /io 167T 2 P 6 
3c 3 hqImoI \ 1/2 / A ^ 1 / 2 


327r 3 P 6 


(") • 


5 > 3.20 x 10 15 (PP) ' (T) 


■ \ 1/2 


(5.59) 


(5.60) 

(5.61) 

(5.62) 

(5.63) 
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(b) For pulsar PSR J0157+6212 the rotational kinetic energy is 
E rot = 2 tt 2 x 10 38 x P ~ 2 = 3.56 x 10 38 J. 

The rate of loss of rotational kinetic energy is 

E mt = 4vr 2 x 10 38 T = 5 . 7 i x 10 25 W. 

P 6 

The pulsar’s characteristic age is 

T = P = 6.2 x 10 12 s = 2.0 x 10 5 y. 

2 P 

The pulsar’s minimum polar surface magnetic field is 

B min = 3.20 x 10 15 (PP) V2 = 2.14 x 10 9 (T). 
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